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Abstract 

We have previously discussed the classical diffusive system of the hounded one- 
dimensional multitrap using the transfer-matrix method which is generally applied for 
studying the energy spectrum of the unbounded quantum Kronig-Penney multibarrier. 
It was shown, by this method, that for certain values of the relevant parameters the 
bounded multitrap array have unity transmission and a double-peak phase transitional 
behaviour. We discuss in this work, using the same transfer matrix method, the energy 
related to the diffusion through the unbounded one- dimensional multitrap and find that it 
may be expressed in two entirely different ways with different results and consequences. 
Also, it is shown that, unlike the barriers in the Kronig-Penney case, the energies at 
one face of the imperfect trap greatly differ from the energies at the other face of the 
same trap. 

keywords: Kronig-Penney system, Transfer matrix, Imperfect trap 
Pacs numbers: 71.15.Ap, 66.30.-h, 02.10.Yn 



1 Introduction 



The remarkable similarity [U [2] between the Schroedinger and the classical diffusion equa- 
tions have attracted many authors to discuss diffusion limited reactions using quantum meth- 
ods and terminology (see annotated bibliography in [1]). For example, the same methods 



and teminology of transfer matrices [31 HI [5] , which are apphed [31 Hj for discussing quantum 
multibarrier potentials, have been used [6l [7] for discussing the one- dimensional bounded 
imperfect multitrap system [HJ [9l [101 [U] through which classical particles diffuse. 

The imperfect trap, which was introduced in [12] and further discussed by others [T3 | [T ^ [T5]. 
may serve as a model for many physical situations. For example, one may find applications of 
it to rotational diffusion in chemical reactions [16j or to proteins with active sites deep inside 
the protein matrix [IT] or to infinite lattice traversed by a random walker in the presence of 
an imperfect trap [13]. We note that the discussion of the bounded one-dimensional multi- 
trap systems have resulted, for certain values of its parameters [6], [7] , in finding somewhat 
unconventional results. Among these one may count a unity transmission of the density of 
the diffusing particles through the multitrap array [6] or the double-peak phase transition 
recently found [7J in such systems. 

An important aspect of the classical one- dimensional multitrap system which was not fully 
discussed thus far is when its length tends to infinity. The analogous quantum infinite 
multibarrier, which is the Kronig-Penney system [31 SI HE] , have been extensively discussed 
in the literature by many authors and it is known by its famous band-gap energy spectrum 
[31 [T8] which is widely applied in electronics, semiconductors and solid state physics [19j. 

We discuss here the problem of a very large (infinite) one- dimensional multitrap system 
using the same transfer matrix method which were applied for studying the Kronig-Penney 
multibarrier potential [31111 HI]- We, especially, discuss the energy of the diffusing particles 
and apply similar methods as those used for studying the enegy spectrum of the quantum 
Kronig-Penney multibarrier [T8l [T9]. 

By using the transfer matrix method for the unbounded classical multitrap system we ob- 
tain a quadratic characteristic equation the two solutions of which give rise to two possible 
expressions for the energy of the diffusing particles. Each of these two energies has a part 
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which is associated with the left hand face of the trap and another, differently expressed, 
part related to the right hand face of it. The different expressions of each of the two energies 
at the left and right hand sides of the trap causes these energies to greatly differ in value 
at these faces. That is, we show that by merely diffusing through the trap the particles's 
energy enormously changes. All the analytical results are graphically corroborated. 

We note that the energies of the bounded one-dimensional multitrap system were found 
in [YJ to have phase-transitional characteristics for the case in which an external field was 
appended to the system. 

In Section II we apply the transfer matrix method for introducing and discussing the un- 
bounded one-dimensional multitrap system as done in [6l [7] . Note that by using the transfer 
matrix formalism we also use its terminology which usually refers to an A^-array system [1] 
rather than to an infinite array one. We remove this finiteness by letting the number of 
traps N and the total length L of the multitrap to become very large. In the numerical 
part we assign to and L the values of 15000 and 20000 (note that in [4] a finite multi- 
barrier potential composed of a few hundred barriers was used as a model for the infinite 
Kronig- Penney system). In Section III we discuss the energy associated with the diffusing 
particles and find, using Appendices A and B, the appropriate expressions for it. In Section 
IV we calculate the energy for some specific values of its parameters. In Section V we use the 
analytical results of Sections III-IV and those of the Appendices A-B for graphically showing 
the energies as functions of its variables. We show that these variables have certain values at 
which the corresponding energy becomes disallowed such as, for example, when it tends to 
become negative or to assume very much large positive values. Some analytical expressions 
and derivations are shown in Appendices A-B. We conclude with a brief summary. 
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2 Application of the transfer matrix method for the 
unbounded one-dimensional multitrap system 



The one- dimensional imperfect multitrap system is assumed to be arranged along the whole 
positive X axis and the diffusing particles which pass through it are supposed to come from 
the negative side of it. We denote, as in [6l[7j, the total width of all the traps and the total 
interval among them by a and h respectively where a and h tend to become very much large. 
The ratio of 6 to a and the total length a + h oi the system are denoted respectively by c and 
L. As in [SI El we may express a and 6 by c and L as a = (Y^j ^ = '^^^ period of 

the multibarrier system which is is denoted by p. We assume that the multitrap system 
begins at the point x = = j^^- 

The initial and boundary value problem [20] which is appropriate for describing the diffusion 
through the imperfect barriers is [6l [7| 



1) pt{x,t) = Dp^^{x,t), t>0, 0<x<{a + b) 



2) p{x,0) = po + f{x), 0<x<{a + b) 
1 dp{x, 
k dx 



1 dp(x, t) , 

3) p{x,,t) = - '\' t>0, l<i<2A^, 



where p{x,t), pt{x,t) and Pxx{x,t) denote respectively the density of the diffusion particles, 
its first partial derivative with respect to the time t and its second partial derivative with 
respect to x. The dissusion constant D is supposed to have two different values; Di inside 
the traps and De outside them where -De > -Dj [HI [7] . The value of 0.5^^ is the order of 
magnitude of the diffusion constant at room temperature and atmospheric pressure (P. 337 
in [21]). In the numerical part here we have assigned to De and Di the respective values 
of 0.8^^ and 0.4^^. The second equation of the set (fll) is the initial condition which is 

sec sec ^ ' — 

assumed [HI [7] to depend on x through f{x) and on the constant term po- The third equation 
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of the set ([T]) is the boundary value condition at the location of the traps where each trap 
has a finite width. That is, any trap is characterized by the two points along the x axis 
where its left and right hand faces are located. The constant k is the trapping rate (or the 
imperfection constant) which characterizes the degree of imperfection of the traps where the 
ideal trap condition is obtained when k ^ oo. The set dl]) may be decomposed into two 
separate problems as follows [51 E] 



1) Pt{x,t) = Dp^^, t>0, 0<x<{a + b) 



2) p(x,0) = po, 0<x<{a + b) 
1 dp{x, t) 



(2) 



3) p{xi,t) = - 



k dx 



t>0, l<i<2N 



1) pt(x, t) = Dpxx{x, t), f > 0, < x < (a + 6) 

2) p(x,0) = /(x), 0<x<(a + 6) (3) 

3) p(xi,t)=0, t>0, 1 < i < 2A^ 



The sets ([2D and ^ respectively represent the diffusion through N imperfect and N ideal 
traps as may be realized from the third equations of these sets. Following [20] one may write 
the general solution of the set ([I]) as [HI [7| 

p{x,t) = Ap^{x,t) + Bp2{x,t), (4) 

where pi{x,t) and p2{x,t) are respectively the solutions of the problems ([2]) and ([3]). Using 
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the method of separating variables [20] one may find the ideal trap solution [6l [7] as 



P2{x,t) = sm( — )e , 



l<i<2N 



(5) 



The solution pi{x,t) of the imperfect trap problem is given by [H 



pi{x,t) = po er/( 



1< i < 2A^ 



2VDt 



2VDt 



(6) 



where the erf{x) and erfc{x) are respectively the error and complementary error functions 
given by erf{x) = e^'^^du and erfc{x) = 1 — erf{x) = f^e~'^^du. The Xi denote the 
2N faces of the traps. Using the transfer matrix method, as done in [3l H] with respect 
to the Kronig- Penney potential and in [6l [7] with regard to the bounded multitrap system, 
one may write the following equation which relates the two faces of the j-th trap 



V 



A. 



\ 



2j+l 



B- 



2j+l 



J 



•~ left s right \ 
,Xj ) 



11 (x^- ,x 



s left s right \ 
Xj ) 



Ti / •■ left •■ right \ 
12 (X/ ,Xj'' J 

TI / •. left •. rinht\ 
22[Xj J 



A2{j-l)+l 

y 52(j_i)+i J 



1<3<N (7) 



A2j+i and i?2j+i are respectively the imperfect and ideal trap coefficients respectively of the 
j-th trap and y42(j_i)+i and i?20-i)+i are those of the (j — 1) trap. The two-dimensional 



left 



matrix T^^^ at the right hand side of Eq ([7]) relates the left hand face of the j-th trap at x^* 
to its right hand face at x"^^* where x"^^* > Xj^* ■ The matrix elements Tn, ri2, T21 and 
T22 are derived in details in [6l [7] and are given in Appendix A. 



pc 



For a one-dimensional A^ trap system, which begins at the point x = — -jjj^ 
period p one obtains the general transfer matrix equation [6l [7] 



and has a 



^ A ^ 

A2N+1 



2N+1 



^^),pAr)T(^-i)(p(Ar - |^),P(iV - 1)), 
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.„r-(p(i+(^),2p)r™(^,p) 



Each two-dimensional matrix at the right hand side of the last equation is denoted in its 
parentheses by the locations of the left and right hand faces of its corresponding trap. 
Thus, one may realize, for example, that for an array which begins, as remarked, at the 
point X = jj^ = the locations of the left and hand side faces of the A^-th trap are 

2^!v'^* = ~ (T+^) x^^^* = pN and those of the first trap are x''^^^ = = -^^^ and 
^^ght _ _ Note that, as remarked in [6], [7j, all the two-dimensinal matrices at the 
right hand side of Eq ([8]) have the same values for D, t, L and c and differ by only the values 
of X along the positive spatial axis. Performing the N products at the right hand side of Eq 
([8]) one may obtains an overall two-dimensional matrix, denoted T/v, whose elements T/v^, 
7/Vi2, 7^f2l T/V22 may be recursively expressed by 



1 1 
^TVn = - , x )^ Np) = ...= T\ Tii{p{j - ),jp) 

(1 + c) fJi (1 + c) 

'^12 = ^Af-l)i2 = • • • = %i2 = T112 =Tl2 = Q (9) 

Tn,, = V-i)2i T22 (piN - ),Np) + r^N-i),, T2, (piN - , Np) 



1 1 
T;v22 = V-i)22^22b(A^ - ^ ^ x ), Np) = ...= Yl T22{p{i - ^ x ),J» 



Note that whereas T/Vn and T/V22 are each a one-term expression which is constructed from 
products the element 77V21 is an A^-term expression each of them is composed of products. 
Now, using Eq {A2) in Appendix A (see also the second of Eqs Q) one may calculate the 
trace Tr and the determinant Det of the two-dimensional matrix T^^^ at the right hand side 
of Eq dZD 
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+T22{P{] 



:i + c) 



(10) 



+ 



,P3, 



■T22{p{] - 



1 + c 



Using Eqs (fTU]) , and following the analogous Kronig- Penney case [21 HI (TS] , one may write 
the following quadratic characteristic equation of T^^"^ 



1 + c 



+T22{pU - 

■T22{p{j 



1 . 



),Pj)) + Tii{p{j - 
pj) = 



(1 + c)- 



<PJ) 



(11) 



The two roots y^^^ and of the last equation which are the required eigenvalues of T*--^-' are 



[l + c) 



),pj), y 



T22{p{j 



(1 + c) 



),pi) 



(12) 



Now, if the two roots y^i\ 1 < j < N are different as for the case here (see Eqs 

(Ai) and (^44) in Appendix A), the two eigenvectors which correspond to them are linearly 
independent and we may identify, as for the corresponding quantum Kronig-Penney system 

from Eq ([H]) with the following two eigenvectors 



[5], the initial values 



7^(1) 



7^(1) 



At 



,(1) 



.(1) 



(At] 

(aA 



Til 
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[l + c) 



,p) 



T22{ ,^ . , ,P) 



(1 + c) 



(13) 



where T^^^ at the left hand sides of Eqs (fT3!l is the two-dimensional matrix from the right 



hand side of Eq ([7]) for j = 1 and and p at the right hand sides of Eqs (fT3|) are, as 
mentioned, the respective locations of the left and right hand sides of the first trap. For 

( A,,, ^ 

these y±.^ one may identify, as for the Kronig- Penney case [3|,|1], the coefficients 



in Eq ([8]) with the two eigenvectors 



V 



B. 



2N+1 



I 



R+ 

I 



A. 



2N+1 



B. 



2N+1 



(14) 



From Eqs (Ai), (A4), (^45) and {Aj) in Appendix A and from realizing that the variables Xi 



assume either the value of x^^-^* or x"^"^ (see, for example, the following discussion before Eq 
( |3T|) ) one may see that the quantities 1 < j < N are identical and satisfy y^^^ — "'^"^^ 



. right 



y+ 



W rjj-^g other quantities y^l\ 



y+ 



1 < i ^ A^ can be seen to slightly differ from each 



other and one may approximately write y 



(1) 



(2) 

y- 



(N) 

yl '. 



Considering the limit of an infinite multitrap array which is arranged along the whole positive 
X axis we should demand, as for the Kronig- Penney case ^|4j|, that as the number of barriers 
A^ tend to 00 the right hand sides of Eqs (iMl) should not diverge. That is, we require 



\y 



221 



+ 
pc 



p) 



< 1 

< 1 



(15) 



Substituting in the last inequalities for Tn and T22 from Eqs (Ai) and (^4) of Appendix A 
one obtains 

a{D^,^,t)a{Di,p,t) 



a{Di,j^yt)a{D,,p,t) 



< 1 



(16) 
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where a and i] are given respectively by Eqs (As) and (A7) in Appendix A. 

3 The energy of the diffusing particles in the one-dimensional 
multitrap system 

In order to be able to reduce the inequalities (|T6ll- (fT711 to calculable expressions we express 
the parameters a and 77, which were given by Eqs {A^) and {Aj) in Appendix A, in terms 
of the energy E of the diffusing particles. We use for that matter the relevant expressions 
of the energy which were fully derived and discussed in [7] for the multitrap system. Thus, 
using Eqs dl])-©, we can write the energy E as 




(18) 




trap coefficients A{x, D) and B{x, D) are numerically found for the 2N faces of the traps 



X = Xj, j = 1,2, . . . 2N [HI E]. That is, for each j-th trap, one may find, using the transfer 
matrix method, the four pairs (1) A{x^^* ,Di),B{x^^^ ,Di), (2) A{x"^^^ , Di), B{x"^^^ , Di), 
(3) A{xf^\ De),B{xY\ De), (4) A{x: 'Y^\ De),B{xf'^\ D^). The first pair denotes the ideal 
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and imperfect trap coefficients inside the j-th trap at its left hand face. The second pair 
denotes these coefficients inside the j-th trap at its right hand face. The third and fourth 
pairs denote these coefficients outside the j-th trap at its left and right hand faces. Note 
that these coefficients, as well as the variables x^'^^^ and xj^^^, are not independent of each 
other. First, one may realize (see the discussion after Eq ([8])) that x^^^^ and x"^''* are given 



by x^^-^* = p{j — (Yqr^), x"^^* = pj, 1 < j < N so that they are related by 



^left^^right_P^^ 1 < j < iV (19) 

■' (1 + c) 

Second, the transfer matrix method relates the former coefficients of the j-th trap among 
themselves and also with those of the {j + l)-st trap as [6j 



A{xY\ A) = A{h^'^, A), A(x7'^\ A) = ^(X(f+\), A) (20) 
B{xY\ A) = S(a;7^'*, A), 5(^7^'*, A) = S(xjf^\), A) 



As one may realize the from Eq (ITS!) does not have to coincide with Xj. That is, although 
for the same j-th trap each of Xj and Xj denote its two faces, Xj may, for a specific context, 
refers to its left hand face in which case it is written as x^^* whereas Xj may refers in this 
context to its right hand face and is written as x^^^*- 

Now, analogously to the Kronig-Penney case [3[|11[T8], one may turn the inequalities at the 
right hand sides of (ITB]) - (IT7|) to equalities. In such case the left hand sides of (ITB1) - (IT7I) are 
equated to cos(/tp) where k is a real parameter and p is the period of the multitrap which 
is p = jj. In accordance with the analogous procedure of the Kronig-Penney case [3], HI [IE] 
the two eigenvalues of the characteristic equation are related to the same parameter. Note 
that even if one relates the two eigenvalues to different parameters he will obtain the same 
following expressions (PT]) - (151|1 for the energies each of which depends on only one parameter. 
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We use in the following the transfer matrix principal property in which the density (and 
its derivative) at the two sides of any face of each trap are equal [3l HI [6]. This may be 
expressed, for example, for the left hand face of the j-th trap as 



p{De, Xi, t) = A{xf^\ De)a{De, xf\ Xi,t) + B{xf^\ De)p2(^e, xf-''^ Xi, t) = 
= p{Di, xY\ Xi, t) = A{xY\ Di)a{Di, xY\ Xi, t) + B{xY\ Di)p2{Di, xf*, Xi, t) (21) 

Substituting from Eq ( fTSj) for the a's in ( fT6l) and from Eqs ([5]), ( fT8|) and (Aj) for the ij^s in 
( fTTl) one obtains after equating the left hand sides of (fT6 l) -( fT7l) to cos^np) (see the discussion 
after Eq (|20D 1 

E(De,j^,i,,t)t ^ B{-^,D^)p2{D^,jf^yX„t) \ / E{D„p,x„t)t _ B{p,D,)p2{D,,p,x„t) 

E{D^,-^,i:^,t)t B(p£^,D,)p2(A.(i^A,t) \ / E{De,pM,t)t B{p,De)p2{De,pM,ty 
^((ifsy.A)A Mj^,D,) J ■ V A(p,D,)D, A(p,D,) 

cos{kp) (22) 



/ E{D,,-^,x„t)t _ A{jJf^,De)a{De,jf^,Xr,t) \ / E{D,,p,x„t)t _ A{p,D,)a{D„p,x„t) \ 
\ B{jJf^,D,)D, BijJ^,De) ) • B{p,D,) ) _ 

( E{Di,j^,Xi,t)t ^ A{j^,Di)a{Di,j^,Xi,t) \ / E{D^,p,Xi,t)t _ A(p,De)a(De,V,Xt-.t) \ 
\ B{^,D,)D, BijJ^,D,) ) ■ \ B{p,D,)D, B(p,D,) ) 

= cos{kp) (23) 



The functions p2 and a are given respectively by Eqs and (A^) in Appendix A and use 
is made of the relation P2{D, x, Xj, t) = —^ri{D, Xi, t) sm{^) obtained by comparing Eq ([5]) 
with Eq (Ai) in Appendix A. The sine function and the factor ^ which do not depend on the 
diffusion constants Di and Dg are cancelled in Eq (!23l) . As realized from the last equations 
there are four energies related to the trap; E{De, -p^; Xi, t), E{De,p, Xi, t), E{Di, -p^, Xi, t), 
and E{Di,p, Xi, t). But, as seen, one may reduce the number of the energies related to each 
trap to two since using Eqs f|T8l) and fl2T]) one may obtain the following expressions which 
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relate the energies at the two sides of each trap 



E{D,,xf\x„t)^ = E{D„xf\x,,t)^, 1<J<N (24) 
E{D,,x^'^^\x,,t)^^ = E{D,,x7^\x,,t)^^, 1<3<N 

In the following we use Eqs fl22|) - flM|) for finding the two energies E{De, p, Xi, t) and E{De, -p^: 
which are respectively the energies at the right and left hand faces of the trap. Thus, using 
Eqs (12^ we may rewrite Eqs (p2|) - (l23|) as follows 

(^E{De,jJ^yX,,t)t-DeB{jJ^,Di)p2{Di,jJ^,x,,t))iE(D<,,p,x^^^^ 

cos(kp) (25) 



[E{De,j^M,t)t-DeA{j^,De)a{De,-^,Xi,t))iEiDe,p,Xi,t)t-DeA{p,Di)a{Di,p,Xi,t)^ 

(E{De,j^,x,,t)t-DeA{j^,Di)a{Di,j^,x,,t))-{E{D^,p,Xi,t)t-D^A{p,D^)a{De,p,x 



D 

■ (l+c) ' 
B{p, De) 



■ cos(kp) 



(26) 



The two quadratic equations fl25l) - fl26l) were simultaneously solved in Appendix B for the 
energies E{De,p, Xi, t) and E{De, (]^) and two separate solutions were found for each 
(see Eqs {B^y^Bg) in Appendix B). For E{De,p, Xi,t) we find the two solutions of 



E^{De,p,Xi,t) 



t{XiX4 — X5)y3 — (X1X2X4 — X3X5)Yi 

{tX;-tXiX2)Yi-t^{l-Xi)Ys 



E {De,p,Xi,t) = — 

J-1 



(27) 
(28) 
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And for E{De, (j^j ^i, t) we find the two solutions of 



i?n^e,7Y^,X.,t) = ^ (29) 



E-iD,,-^^^,x,,t)= (30) 

_ {X3-XiX2)({tX3-tXlX2)Y2HXlX2X4-X3X5)Yi)-{l-Xi)ii'{X3^XiX2)(Y4-Ys)+t{XiX2X4-X3Xs)Y3) 

(X3-XiX2)((t2(i-Xi)y2+*(XiX4-X5)Yi)-(i-Xi)(t3(i-Xi){y4-y5)+t2(XiX4-X5)y3) 

where the quantities Xi, X2, X3, X4, XsandFi, Y2, Y3, Y4, Is are given respectively 
by Eqs (-Bi) and (S4) in Appendix B. The energies E~^{De,p, Xi,t) and E^{De,p,Xi,t) from 
Eqs fl271)- fl28l) are for the right hand side of the trap and the energies E~^{De, (j^, Xi, t) and 
E~{D^, Xi, t) from Eqs ([29])-(l30]) are for the left hand side of it. 

In the former expressions of the energies the variable Xi must coincides with either (]^^ 
or p. Thus, when Xj = one have to discard the solutions £'^(De, -p^, -p^, ^) and 

E^{De,p, ji^i since in this case one obtains from Eq ([5]) and from Eqs (Bi) in Appendix 
B X2 = = 0. In this case the energy E^{De, (j^) (i^i ^) from Eq fl30|) vanishes whereas 
the energy E^ {De,p, j^^j'^) from Eq fl28|) remains at the value of ||. This implied the 
unreasonable conclusion that the passing particles have no energies at the left hand face of 
the trap before they diffuse through it whereas at the right hand face of it they have non- 
vanishing unaccountable energies. Thus, for Xi = only the energies E^{De,p, j-^^,t) 
and E~^(De, j^^^ life)' rnust be considered and they are given by 

^ ^"^-^'(TT^'^^^ t(X,-l) ^''^ 

The second case is that of Xi = p for which one obtains from Eq ([5]) and from Eqs [Bi) in 
Appendix B X4 = X5 = 0. In this case the energy E^{De,p,p,t) from Eq (|27|1 vanishes 
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whereas the energy E^{De, ^^^^y P^ t) from Eq (129!) remains in the value of This also 
could not be accepted since it implied that the particles suddenly stop diffusing after the 
first trap whereas we are concerned with the diffusion along the entire multitrap system. 
Thus, the energies E~^{De,p,p,t) and E~^{De, -^Y^,p,t) must be discarded and we have to 
take into account only the energies E~{De,p,p,t) and E^{De, p^iP, which are given by 

E-{D,,p,p,t) = ^ (34) 

If c becomes very large so that c >> 1 one may realize from Eqs (Bi) in Appendix B and Eqs 
(E]), (EH) and (133|) that the energies E^{De,p, j^^,t) and E^{De, j^^^P^i) tend to zero. 



4 Calculation of the energies (1311)- (1341) for specific val- 
ues of H.p. 

The expressions for the energies E+{D,,p, ^,t), E+iD„ ^,t), E' (D^, ^,p, 

E^ {De,p,p,t) should now be evaluated as functions of up as done for the quantum Kronig- 
Penney system [3l IH |18]. But before doing that we show that the expressions fl3Tl) - flMl) 
become simplified for certain values of up. Thus, for /tp = | + nir, n = 0, 1,2, . . . one 
have cos(kp) = and from Eqs {B4) and the first of Eqs {Bi) in Appendix B we have 
Yr, = Xi = and also the second terms of Yi, Y2 and Y^ vanish as well. Thus, the energies 
(EU-dMl) become 

PC X, _ DeB{p,D.)p2{D.,P,jftyt) 
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(cos(Kp)=0)V-e, ^^^^y ^^^^y 

■ eicii{k'^Det)erfc{k^/D~t) 



t 



t 

(36) 




pc 



(l+c) 



pc 



p,t) 



(37) 



(l + c) 



t t 



E, 



(cos(kp)=0) 



{De,p,p,t) 




■ exp{k^Dit)erfc{k>/D^) 



(38) 



For obtaining Eqs 0351) and (13 7p we respectively use the fifth and third of Eqs (-Bi) in 
Appendix B and for Eqs (136!) and (!38l) we respectively use the third and second of Eqs (S4) 
in Appendix B. Use is also made of Eq (A^) in Appendix A and the first of Eqs (-B4) in 
Appendix B. 

Another kind of points which draw special attention is Kp = arcjcosi ^/ ^^pc'^ '^p t)) for which 

(l + c) ' 

one obtains from the first of Eqs (Bi) in Appendix B Xi = 1. At these points the energies 
E^{De,p, -^j^,t) and E^ (D^, -^j^,p,t) from Eqs fl3Tl) and (!33|) must be discarded since 
they tend to infinity and this can not be accepted on physical grounds. Note that using 
the transfer matrix method one may conclude that the ideal traps coefficients always sat- 
isfy B{x"^^^ , De) > B{x^^* ,Df.) and so the division ^{ ^^^^ ''d ) greater than unity which 

(l + c) ' 

implies that the quantity Yi, as defined by the first of Eqs (S4) in Appendix B, is always 
positive. This determines, as will be shown, the values and the graphical form of the energies 



(l+c)' (l+c) 



pc pc 



,t) and E'{De,p,p,t) from Eqs (|32|) and (IMl). 
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5 The energies as functions of Kp, c, k, and t. 



As seen from Eqs (Bi) and (S4) in Appendix B, (^45) in Appendix A and from Eqs ([5])- 
the energies fl3T|) - fl34|) critically depend upon the ratio c, the trapping rate k and the 
time t. Also, one may conclude from the analytical form of the expressions fl31l) - fl34l) and 
from the Panels of Figures 1-4 that the energy £'+(De,p, -jj^jt) from Eq (13T]) corresponds 
to E-{De,j^yP,t) from Eq ([33D and E+{De, j^yt) from Eq ^ corresponds to 
{Df,,p,p,t) from Eq flM|) . That is, E^{De,p, j^^,t) and E^ {Df,, j^^,p,t) are expressed 
only by the ideal trap expressions from Eqs (Si) in Appendix B and E~^{De, -p^i (i^i ^) a^id 
E~{De,p,p,t) are given only by the imperfect trap expressions from Eqs {B4) of Appendix 
B. The corresponding energies E~^{De,p, i^^jt) and E~ {D^, -ij^,p,t) as functions of Kp 
are characterized with a behaviour which causes them to abruptly change their values in a 
rather jumpy and discontinuous way (see, for example, the Panels of Figures 1-2). As one 
may assume these abrupt changes in the energies fl3T]) and fl33|) are related to the values 
of Kp for which {Xi) in their denominator is close to 1. Note that the nonzero values of 
these energies may be negative in which case they can not represent real energies since we 
discuss here only kinetic energies as realized from Eq (fTSl) . The second corresponding pair 
of energies E^(De, j^^, -^-^^^t) and E~{De,p,p,t) are characterized as steeply increasing 
with up for very small values of it and at Kp ^ 0.05 they become constant (as functions of 
Kp) for all Kp > 0.05 (see, for example, the Panels of Figures 3 and 4). Also, in contrast to 
the former pair, these energies are always positive. 

From the Panels of Figures 1-4 one may realize that, for the same values of c, k and t, the 
nonzero values of the energy E^{Df,, jj^,p, t) are generally greater by several orders of mag- 
nitude from the other three energies (see, for example. Panels 1-2 of Figure 4 which show a 
giant difference of 10"^^ between E~{D(,, -p^iP, i) and E'^i^De^p, -p^, t)). Thus, although, as 
mentioned, the energies E~ {De,p,p, t) and E'{De, 7^7^, P, t) correspond in analytical expres- 
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sions and graphical form as functions of Kp to the respective energies E'^{D^, Ju^)' (h^' 
and E~^(Df.,p, -p^^, t) they greatly differ in value. This may be seen from Panels 1 and 2 of 
Figure 1 which, respectively, show the energies E~{De, jY^,p,t) and E~^{De,p, -^^^,t) as 
functions of Kp using the same 20 different values of the ratio c for each Panel and the same 
t and k for all the graphs shown. The 20 values of c are c = 0.1 + n • 0.4, n = 0, 1, 2, ... 19 
and the values of t and k for all the graphs shown in the two Panels are t = k = 2. Note that 
although the energies shown in Panels 1 and 2 look similar as functions of up the nonzero 
values of these energies differ by as much as 10^. As mentioned, we should consider only the 
zero or the positive parts of the graphs as representing real kinetic energy. 

In Panels 3 and 4 of Figure 1 we show enlarged views of the small sections in the respective 
Panels 1 and 2 just to the right of the point Kp — 0. Note the similarity between these 
energies even at this small resolution and also note that despite this similarity the nonzero 
parts of the energy in Panel 3 are about 0.5 • 10^ erg whereas those of Panel 4 are about 
0.4 erg. The largest hooked negative graph corresponds to the smallest value of c and as c 
increases the other hooked positive and negative graphs arc added. The larger c becomes in 
Panels 3-4 the corresponding energies become smaller and tend to be densely arrayed around 
zero. This means that the larger is the interval between the traps compared to their width 
the kinetic energy of the diffusing particles tends to decrease to zero. 

The same similarity in graphical form and same large differences in values may be shown 
for the same energies from Figure 1, as functions of Kp, but now for different values of 
the trapping rate k. This is seen in Panels 1 and 2 of Figure 2 which respectively show 
the energies E~{De, j^^^Pi't) and E^{Dg,p, -^j^yt) as functions of Kp using the same 20 
different values of k for each Panel and the same t and c for all the graphs shown. The 20 
values of k for each panel are k — 0.1 + n- 0.4, n 1, 2, . . . 19 and the values of t and c for 
all the graphs in the two Panels are t = 2 and c = 1. As for Panels 1 and 2 of Figure 1 the 
differences between the nonzero parts of these energies amount to about 10^ although they 
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Panels 1 and 2 respectively show the energies E {De, j^^^,p,t) and 



t) as functions of up for the 20 values of the ratio c = 0.1-|-n-0.4, 



n — 



0,1,2,.. .19. Panels 3 and 4 show a hight resolution of the respective neighbour houds in 
Panels 1 and 2 just to the right of the point up = 0. Note that although Panels 1 and 2 as 
well as Panels 3 and 4 are similar in form they greatly differ in the nonzero values of their 
energies. Both the trapping rate k and the time t have the values oi k = t = 2 for all the 
graphs of the four Panels. The units of the energies are in ergs. 
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look similar in external form. Note that actually the energy in Panel 2 tends to zero. Panels 

3 and 4 of Figure 2 respectively show enlarged views of the respective neighbourhouds from 
Panels 1 and 2 about the point np — 25. One may note the similarity between these energies 
even at this small resolution. Also, one may note that despite this apparent similarity the 
nonzero parts of the energy in Panel 3 is about 10^ whereas the corresponding ones in Panel 

4 tend to zero. 

For both Panels 3 and 4 the graphs with the large dense hooked positive parts correspond to 
the smallest values of k which means, as one may assume, that the smaller is the trapping 
rate of the traps the larger is the energy of the diffusing particles. As k increases the 
corresponding graphs become negative and they tend to zero for large enough k. That is, 
the more k grows which means that the larger becomes the trapping rate of the trap the 
more restrained and blocked become the diffusing particles in their passage through it. This 
is demonstrated through the vanishing of the positive allowed parts of the energies for large 
k and their tendency to the zero value. 

In the Panels of Figures 1-2 we compare for different values of c and k the two corresponding 
energies E~{De, -(^^,P, t) and E^{De,p, (]^) t) as functions of up. We now discuss the sec- 
ond pair of corresponding energies E~ {Dg^p,p^t) and E'^ {Dg , , , t) . Compared to the 
enrgies E~{De, jY^-,P-,t) and E^{De,p, (j^^^) from Figures 1-2 the energies E~{Df,,p,p,t) 
and E^{De, j^f^, (i^j ^) always positive and they are generally constant with np. Pan- 
els 1 and 2 respectively show three-dimensional surfaces of the energies E~{Df.,p,p,t) and 
E^{Dc, (1^1 (1^' ^) ^ functions of Kp and c and for the values ol k — t — 2. Note that the 
energy E~^{De, jif^, ' ^) Panel 2 does not depend at all on either np or c and has the 
rather small constant value of 0.088 erg. The energy E~{Df.,p,p,t) at Panel 1 is constant 
for all /cp's and depends only slightly on c as seen from the small depression of the surface 
at small c which causes it to be slightly distorted from the planar form of Panel 2. Panels 3 
and 4 respectively show three-dimensional surfaces of the same energies from Panels 1 and 



20 




Figure 2: Panels 1 and 2 respectively show the energies E {De,j^^,p,t) and 
E'^{De,p, j^^,t) as functions of Kp for the 20 values of the trapping rate k — 0.1 + n • 
0.4, n = 0,1, 2,... 19. Panel 3 and 4 show a hight resolution of the respective neigh- 
bourhouds in Panels 1 and 2 of the point Kp = 25. As is the case for the Panels of Figure 
1 one may note that although Panels 1 and 2 as well as Panels 3 and 4 are very similar in 
form nevertheless they greatly differ in the nonzero values of their energies. The ratio c and 
the time t have the respective values of c = 1 and t = 2 for all the graphs of the four Panels. 
The energies are in units of ergs. 
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Figure 3: Panels 1 and 2 respectively show three-dimensional surfaces of the energies 
E~{De,p,p,t) and E^{De, p^^, jif^^t) as functions of up and c and for /c = t = 2. Note 
that the energy E^{Df,, -(^+c)i (i+c) ' ^) Panel 2 does not depend on either np or c and has the 
constant value of 0.088 erg. The energy E^{De,p,p,t) at Panel 1 depends only slightly on 
c, is constant for all rep's and have a maximum value (for = 2) of E^^^{Dg,p,p, t) — 1 erg. 
Panels 3 and 4 respectively show three-dimensional surfaces of the same energies from Panels 

1 and 2 but now as functions of Kp and k and for c = 1 and t = 2. Note that these energies do 
not depend at all on Kp and vary with k to the maxima (for c = 1) of E^^^{De,p, p, t) = 15 erg 
and £'+„^(De, (j^, (j^, t) = 0.23 erg. Note also that the energies E~{De,p,p,t) and 
E~^{De, (1^, (i^) drop respectively to zero ai k ^ 4.8 and k ^ 6.3. 

2 but now as functions of Kp and k and for the values c = 1 and t = 2. Note that these 
energies, as in Panels 1 and 2, do not depend on K.p and vary with k to the maxima (for 
c = 1) of E-^^{De,p,p,t) = 15 erg and £;+„,(L'e, (i^, (i^,^) = 0.23 erg. Note also that 
as the trapping rate k grows the energies E^{De,p,p,t) and E'^{De, j^^^i (]^'^) decrease 
in value and become zero at the respective values of A; ~ 4.8 and k ~ 6.3. This result is 
expected since as the trapping rate grows the traps become more effective in blocking the 
diffusing particles. 
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As seen from the Panels of Figures 1-3 the diffusing particles's energy considerably changes by 
merely passing through the trap. Thus, referring to the pair E~{De, -p^, i), E^iDe,p, p, t) 
one may realize the large change in kinetic energy the diffusing particles goes through upon 
passing from the left hand side to the right hand side of the trap. For example, comparing 
Panel 1 of Figure 1, which shows the energy E^{De, j^^iPii) at the left hand side of the 
trap to Panel 1 of Figure 3 which shows the energy E^{De,p,p,t) at the right hand side of 
this trap for the same values of c, k and t one may realize that the particles's nonzero values 
of the energy changes upon diffusing through the trap from E ^ 10^ erg to i? ~ 1 erg. 
These large differences may be realized again by comparing Panel 1 in Figure 2, which shows 
the energy E~{Df., jY^,p,t) at the left hand side of the trap to Panel 3 of Figure 3 which 
shows the energy E^{De,p,p,t) at the right hand side of it for the same values of k, c and 
t. As seen, the particles's nonzero values of the energy decreases upon passage of the trap 
from E ^ 10'^ to E ^ 15. Thus, one may conclude that by diffusing through the traps the 
particles lose a huge amount of the energy they possess before the diffusion. 

The time evolutions of the energies from Eqs fl3T]) - flMl) as functions of up reveal in a more 
pronounced way the mentioned large differences in the nonzero values of the energies. This 
is demonstrated in the first two Panels of Figure 4 which show the energies E^{De, j^^,P, t) 
and E'^{De,p, -p^; t) for the 60 different values of t = 1 + ri ■ 0.5, n = 1, 2, ... 59 in each 
panel and for c = 2 and k = 1 for all the graphs shown. Note the giant differences of about 
10"^^ between the nonzero values of the energies in Panels 1-2. Continuing to increase t 
causes the energy in Panel 1 to grow (not shown) even beyond 10^° erg. Since these energies 
are not physically possible we conclude that there exist points along the up in which the 
energies E~{De, jj^,P, t) are not allowed for large values of the time t. In Panel 3 we show 
a three-dimensional surface of the energies E~{De,p,p,t) as function of Kp and the time t. 
Note that it is constant with Kp and decreases to zero, in contrast to the energy from Panel 
1, as t increases. Panel 4 shows the energy E~^{De, jif^^ liTc)^ function of Kp, for the 
20 values of t = 1 + n ■ 0.5, n = 0, 1, ... 19 and for c = 2 and A; = 1 for all the graphs. The 
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dense line just above the abcissa axis denote the higher values of t for which the constant 
values of the energy (as function of Kp) tend, like those of Panel 3, to zero. 



From the discussion thus far one may realize that generally the nonzero values of the energy 
E~{De, j^^^,P,t) are greater by several order of magnitudes from the other three energies 
as shown by comparing the Panels of Figures 1, 2 and 4. These great differences are further 
pronounced for increasing values of the time as shown in the Panel 1 of Figure 3. But that 
is no more so when the time decreases as turns out (not shown) when the energies were 
calculated at small times. Thus, for example, decreasing the time from t ^ 30 to t ^ 0.2 
causes the nonzero values of the energy E^{De, j^^tVi^) to decrease from about 10^'' erg 
(see Panels 1 of Figure 4) to about 100 erg. Likewise, the energy E~{Df.,p,p,t) decreases 
from about 15 erg (see Panel 3 of Figure 3) for t = 2 to 3 erg for t = 0.2 (not shown). 



6 Concluding Remarks 



We have discussed, using the transfer matrix method, the energy of the particles which diffuse 
through the unbounded one-dimensional multitrap system. The classical initial and bound- 
ary value problem related to diffusion through an imperfect trap was adapted to apply to an 
infinite array of similar traps as done in the sets ([I])-([2])- Following the conventional trans- 
fer matrix procedure, which is used for discussing the quantum Kronig-Penney multibarrier 
array, we obtain a similar matrix equation (Eq ([S])) which relates the imperfect traps across 
the whole array. Using, as for the analogous quantum multibarrier system, the periodicity 
of the array we obtain a quadratic characteristic equation (Eq f|TT]) ) for the two-dimensional 
matrix which relates the two faces of the general j — th trap. We solve this equation for 
the involved eigenvalues of this matrix and impose upon them the finitness condition at the 
limit at which the number of barriers becomes very large. As a result two inequalities 
(|T6i)-( fT71) are obtained which are the central expressions from which we derive the appropri- 
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Panels 1 and 2 respectively show the energies E {De, ^^^^^^ 
jj^,t) as functions of up for the 60 values of the time t = 1 + n ■ 0.5, 



p, t) and 



n 



0,1,2, ...59. One may realize that as the time grows the nonzero values of the energy 
E~{Dg, j^^,p, t) steeply increase. The ratio c and the trapping rate k for all the graphs of 
Panels 1-2 are c — 2 and k — 1. Compared to the energy from Panel 1 which grows with t 
the energies E~ {De,p,p, t) and E~^{De, -p^i ^) ^^own respectively in Panel 3 and 4 de- 
crease with time to zero. Panel 4 is drawn for the 20 values of t = H-n-0.5, n = 0, 1, 2, ... 19 
and for c = 2 and k = 1. The upper lines in Panel 4 fit the small values of t and the lower 
lines fit the large values. The energies are given in units of ergs. 
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ate kinetic energies of the diffusing particles. Writing the matrix components Tii(-^j^,p) 
and T22(-p^,p) in these inequahties in terms of the appropriate energies (see Eqs (fT5l) - (!23l) 
and discussion there) and using the properties of the transfer matrix method (see Eqs ( |T8|) . 
fl2T|) and flMl) ) we obtain two simultaneous equations involving two energies. We have found 
that each of these two energies is composed of two parts; one is related to the left hand 
face of the trap and the second to its right hand face. It is also found that the two parts 
of each of the two energies differ greatly from each other not only in value but also in the 
way they are expressed as functions of the related variables c, k, t and Kp. That is, as seen 
from Eqs fl3T|) - (IM|l one part of each of these two energies is expressed in ideal trap terms 
only and the second part in imperfect trap terms only. These differences entail the results 
that by merely diffusing through the trap the particles's energy totally changes as realized 
from the appended figures. Moreover, there exist great variations not only between the two 
parts of the same energy but also between different sections (along the Kp axis) of the same 
part itself. For example, in the respective Panels 1 of Figures 1, 2 and 4 we have found that 
the energy E^{Df,,x^^-^^,p,t) at the left hand face of the trap assumes values which greatly 
varies even in very short ranges of up. 

As discussed in Section IV the points along the up axis at which the energies may assume 
unexpected, or even disallowed values, are related to the following two kinds of points; (1) 
= I + nn, n = 0, 1, 2, . . . (2) up = arc(cos( /ie/t ' \ ))- Another important variable 
which entails a large changes in the values of the energies is the time t as realized from 
the Panels of the appended Figures (see, especially. Panels 1-4 of Figure 4). The analytical 
expressions obtained are corroborated by the different attached figures. 

As noted, the analogous discussion of the quantum Kronig-Penney multibarrier entails the 
finding of points along the corresponding Kp axis at which the energy is disallowed. Here, 
for the classical imperfect multitrap we have found corresponding disallowed energies which 
take the form of either negative values for the kinetic energy or of a discontinuous change of 
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this energy from zero to enormous positive values as in Panel 1 of Figure 4. The quantum 
band-gap structure found in the Kronig-Penney multibarrier array have been turned out to 
have great applications in wide areas of solid state physics such as semiconductor devices 
and computer chips. The striking similarity in the forms of the Schroedinger and diffusion 
equations as well as the common possibility to investigate and discuss them by the transfer 
matrix method may entail in the future similar successful development for the classical 
diffusive systems. 

APPENDIX A 



A The matrix elements from Eq (^!) 

rpi J- • 1 i. rr /•■left riqht\ rji /--left riqht\ rri /•-left riqht\ j rji /--left riqht\ 

ihe matrix elements lii[Xj i^j )■> -^i^yXj , a;^- ), -t 2i(a^j ); and ^22(2;^ , a;^- ) 

of the two-dimensional matrix T^-'^ from Eq ([7]) are fully discussed and derived in [6l [7] and 
are given by the following expressions 



rr /-left .right\ Oi{Df,,Xj , Xi,t)a{Di, Xj ^ ,Xi,t) . ^ ■ ^ ,r /A^ 
Tii[Xj , Xj- ) = Zl^fr~s . right ' 1 ^ J ^ ^ (^1) 



Tu{x;^\ x;'""') = 0, 1 < J < TV (A2 

rp /-left -right\ _ /Vi^ijXj^ , t) XiD^, X j , Xi , t) 
l2l[Xj ,Xj ) — Pol ^ 



r]{De,x-^ ,t) r]{Di,x-^ ,t) 



oc^yD^^ Xj , Xi^ t^^iylDi^ Xj , Xi^ 
a{Di, xY\ Xi, t)r]{Di, xY\ t) 



)) + 



( / \ riaht ,\ (^3 



a{Di,x-\Xi,t) ri{De,x-^ ,t) 

), l<J<iV 



oiiyDi^Xj J x^^t^^i^D^^ Xj ^X'l^t^ 



a{De,x"^ ,Xi,t)?7(De,x7^ ,t) 
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The parameters a, ^, and r] are given by (we write these expression for D^, and x = x 



■j > 



.left\ 



a{De,x^^ ,Xi,t) = erf{ 




) + exp{k'^Det + k{x- 



left 



Xi))- 



■erfc{k^D~t + 




(A) 




xf\ x„ t)=k expik-'D.t + k{xY' - x,))erfc{k^t + ' {A,) 



Note that in [6l [7] the variables Xi are not subtracted from the variables Xj in the functions 
a and ^. This is because the presence or absence of this subtraction do not affect at all the 
values of the matrix elements Tn and T22 as may be realized from their definitions in Eqs 
{Ai) and (^44) in this Appendix. Also, in P [7] we discuss the whole array of the bounded 
dense multitrap in which case the variables Xi and Xj do not, necessarily, refer to the same 
trap and so this subraction is ignored there. Here, on the other hand, the variables Xi 
and Xj refer to the same trap which represents the unbounded multitrap system and so the 
expression {xj — Xi) should not be approximated to Xj. 



r]{De,Xi,t) 
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APPENDIX B 



B The solutions of the simultaneous Eqs (1251)- (f26l) 



We solve in this Appendix the two Eqs fl25|) - fl26|) for the energies E{De, j^^,Xi,t) and 
E{De^p,Xi,t). We begin by solving Eq ( 1251) for E{De, j^^^,Xi,t) in terms of E{De,p,Xi,t) 
and then we substitute this solution in Eq (l26l) and solve it for E{De,p, Xi,t). In order not 
to be involved with cumbersome expressions we define the following quantities 

= Aip,D,) • ^"^('^^^ 

Xs = DeB{-^,De)p2{D,,-^,Xi,t) (Bi) 
(1 + c) (1 + c) 

X4 = DeB{p, De)p2iD^,p,Xi,t) 

X5 = DeB{p, Di)p2{Di,p,Xi,t) 

Substituting the last quantities in Eq (1251) we obtain 

e{l - Xi)E{D,, X,, t)E{D,,p, X,, t) - t(E{D,, x,, t){X, - Xai) + 

(1 + c) V (1 + c) 

+E{De,p, Xi, t){Xs - + X3X5 - X2X4X1 = (S2) 

Solving the last equation for E{De, j^^: ^ii t) we obtain 

pc . ^ X.X^jX^ - tE{D„ p, x„ t)) - X^jX, - tE{D,,p, x,, t)) 
[ + tXiiX,-tEiD,,p,x^,t))-t{X,~tEiDe,p,Xi,t)) ^ 
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We may now substitute the last expression for E{De, j^^^Xi^t) in Eq (126!) and solve it for 
E{De,p, Xi,t). But before proceeding we define the following quantities 



,2, B{p,D, 



Y2 = tDe [ j ^^lc Di)ct{Di,p, Xj, t) - A{j), De)a{De,P; Xj, t) cos(/tp) 

^B{j^,De) (1 + c) (1 + c) (1 + c) 

■a{Di,-— — -,Xi,t)cos{Kp)) {Bi 
(1 + c) / 

Y, = Dl f^^:^'l A{p, DMD^,P, t)A{7rr-y DMDe, 7-^, t) 
B{jf^),De) (1 + c) (1 + c) 

Y5 = DlA{p, De)a{De,p, Xi, t)A{ Di)a{Di, J^^ Xi, t) cos(kp) 

(1 + c) (1 + c) 



Substituting from Eqs {B^)-{B4) in Eq (1261) and rearranging we obtain the following quadratic 
equation for E{De,p,Xi,t) 

E\D,,p, t) {{tX^ - tXiX2)Yi - t\l - Xi)F3) + E{D,, p, x„ t)- 

■[{t\l - X,){Y, - Y,) - {tX^X^ - tX,)Y, - t(X3 - X,X2)Y, + {X,X^X2- 

-XsX,)Y,^ - (X4X1X2 - ^3X5)^2 + tiX.Xi - X,){Y, - F5) = (55) 
The two solutions of the last quadratic equation are 

p+m -n ^ _ t(XiX4 - X5)y3 - (X1X2X4 - X3X5)yi 

^ (^-^'^-^) - itX,-tX,X,)Y,-t^{l-X,)Y, ^^'^ 
E-(D f)- (^3-XiX2)y2-t(l-X0(F4-r5) _Y, 

^""-^'"^''^ - (X3-x,x,)y,-t(i-x,)y3 " ^""'^ 
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where the last result for E~{De,p^ Xi, t) is obtained by using Eqs (-B4). The two expressions 
from (Bq)-{Bi) are the energies at the right hand side of the trap. The corresponding 
energies E^{De, j^^,Xi,t) at the left hand side of it may be obtained from Eq (B3) by 
substituting in it for E^{Df.,p, Xi, t) from Eqs {BQ)-{eB'j). Thus, using Eq (Sg) one may find 

(1^, Xi, t) as 

E+{D^,-^ Xi,t)^^ (Bs) 
(1 + c) Yi 

where use is made of the first and third of Eqs {B4). The second energy E~{De, ■^^^,Xi,t) 
is obtained by substituting from Eq (S7) in Eq (-B3). 



(X3-XlX2)((tX3-tXiX2)y2+(XiX2X4-X3X5)yi)-(l-Xi)(t2(X3-XiX2)(y4-y5)+t(XiX2X4-X3X5)l3) 

(X3-XiX2)((t2(i-Xi)y2+t(XiX4-X5)yi)-(i-Xi)(ta(i-Xi)(y4-y5)+t^(^i^4-X5)y3) 
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